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ABSTRACT 



Context. PSR J1713+0747 is a binary system comprising millisecond radio pulsar with a spin period of 4.57 ms, and a low-mass 
white dwarf (WD) companion orbiting the pulsar with a period of 67.8 days. Using the general relativistic Shapiro delay, the masses 
of the WD and pulsar components were previously found to be 0.28 ± 0.03 M Q and 1.3 ± 0.2M o (68% confidence), respectively. 
Aims. Standard binary evolution theory suggests that PSR J1713+0747 evolved from a low-mass X-ray binary (LMXB). Here, we 
test this hypothesis. 

Methods. We used a binary evolution code and a WD evolution code to calculate evolutionary sequences of LMXB s that could result 
in binary millisecond radio pulsars such as PSR J1713+0747. 

Results. During the mass exchange, the mass transfer is nonconservative. Because of the thermal and viscous instabilities develop- 
ing in the accretion disk, the neutron star accretes only a small part of the incoming material. We find that the progenitor of PSR 
J1713+0747 can be modelled as an LMXB including a donor star with mass 1.3 - 1.6M Q and an initial orbital period ranging from 
2.40 to 4.15 days. If the cooling timescale of the WD is 8 Gyr, its present effective temperature is between 3870 and 4120 K, slightly 
higher than the observed value. We estimate a surface gravity of Log(g) at 7.38 - 7.40. 

Key words, binaries: general - stars: evolution - stars: white dwarfs - pulsars: individual (PSR J1713+0747) 



1. Introduction 

Since the first X-ray source Sco X-l was discovered 

dGiacconi et al. 19621 , it was proposed that the strong 

Galactic X-ray sources were accreting neutron stars or 



black holes in binary systems ( 


Zeldovich & Guseynov 1966 


Novikov & Zel'Dovitch 1966, Sh 


klovskii 1967). These sources 



have significant higher flux in X-ray than other wavelengths. 
Over 90% of strong Galactic X-ray sources can be classified in 
two kinds as follows: high-mass X-ray binaries and low-mass 
X-ray binaries (LMXBs). The latter consists of a neutron star or 
a black hole and a low-mass donor star (S; 1 - 2M ). 

The mechanisms that drive mass transfer in LMXBs de- 
pend on the initial separations of the binary components 
(Bhattacharya & van den Heuvel 1991). In narrow systems with 
an orbital period (P rb) less than two days, mass transfer is 
driven by angular momentum loss due to gravitational radia- 
tion and magnetic braking. However, mass transfer is driven 
by the evolution of the secondary in relatively wide LMXBs 
(P or b ^ 2 d). In LMXBs, mass accretion onto neutron star 
causes a magnetic field decay and spins it up to a short spin 
period (;£ 20 ms). It is believed that circular binaries in- 
cluding millisecond radio pulsar and He or CO white dwarfs 
(WDs) are formed as a result of long-term mass accretion 
flAlparetal. 1982| IStairs 20041 . The discovery of the first mil- 
lisecond X-ray pul sar in LMXB SAX 1808.4-365 8 strongly sup- 
ports this scenario ( Wijnards & van der Klis 1998| l. 



Through a survey for millisecond pulsars with the 305 m 
Arecibo radio telescope, poster et al. (1993)| have discovered 
PSR J 17 13+0747. This pulsar has a spin period of 4.57 millisec- 
ond. It is in a circular orbit of 67.8 days and has a low-mass WD 
as companion. It has a weak magnetic field of 1.9 x 10 8 G and a 
spin-down age of 8.9 +1.0 Gyr, so is an older neutron star. By 
detecting the general relativistic Shapiro delay, the neutron star 
and WD were set up as the low limits for the masses of 1 .2M 
and O.27M , respectively dCamilo et al. 1994K 

Through observations for 12 years between April 1992 and 
May 2004, 343 pulse times of arrival for PSR J17 13+0747 were 
recorded with uncertainties down 200 nanosecond. These timing 
data can result in the measurement of the Shapiro delay, which 
can constrain the mass of the WD Mo and the inclination of the 
orbit i. The observational data constrained the masses of the pul- 
sar to 1.3+O.2M an d the WD to 0.28+0.03M o (68% confidence, 
|Splaver et al. 2005| ). 

It is clear that PSR J1713+0747 is a recycled pulsar, be- 
cause it gained mass and angular momentum from its donor 
star, the progenitor of the WD (Alparetal. 1982 1. Assuming 
the mass of the WD is equal to the helium core mass of 
its progenitor, binary evolution theory has presented a pe- 
culiar relation between the orbital period P 01 -b and the mass 
of 



Tauris & Savonije 1999 



WD M 2 flRappaport et al. 1995] 
Podsiadlowski et al. 2002l l. According to this relation, Tauris & 
Savonije (1999) obtained the range of the WD mass for PSR 
J1713+0747, O.31M < M 2 < O.34M . Later on, M 2 was 
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constrained in the range of 0.30 - O.35M dPfahl et al. 2 002 ). 
After adopting the predicted value by the theoretical relation 
of P or b - M2, the pulsar has a somewhat higher mass, Mi = 



1.53^M (68% confidence) (Spl aver etal. 2 005). 



Assuming that the neutron star accretes 10 % of the mate- 
rial coming from the donor star, Ben venuto et al. (2006)| have 
calculated the evolution of PSR J17 13+0747 progenitor to fit 
its component masses and the orbital period. Their results show 
that the progenitor of PSR J 17 13+0747 may be an LMXB with 
initial masses of Mi = 1.4M and M2 = 1.5M Q and an initial 
orbital period off; = 3.10 d (Z=0.02) or P; = 3.05 d (Z=0.01). 
Adopting a semi-analytic approach, Chen et al. (2006) have also 
performed binary evolution calculations to account for the com- 
ponent masses and the orbital period of PSR J1713+0747. 

In this paper, where we employ a detailed stellar evolution- 
ary code and a WD evolutionary code, we attempt to simulate 
the evolutionary history of the progenitor of binary millisecond 
radio pulsar PSR J 17 13+0747. The structure of this paper is or- 
ganized as follows. In section 2, we describe the stellar evolu- 
tionary codes, input physics in the binary evolution calculations, 
and the WD evolutionary code. Numerically calculated results 
for the evolutionary sequences of binary system resulting in the 
birth of PSR J 17 13+0747 are presented in section 3. Finally, we 
discuss and summarize our results in section 4. 



2. Stellar evolutionary code and input physics 

By employing an updated version of the stellar evolution- 
ary code develo ped by |Eggleton (197T)| |Eggleton (1972)| (see 
IHan et al. 19941 IPols et al. 19951 ), we investigated the evo- 
lution of LMXBs. The LMXBs consist of a main se- 
quence donor star with a mass of M2 and a neutron star 
with a mass of Mi. Radiative opacities are those from 
and 



Rogers & Iglesias (1992) and at low-temperature molecular 
opacities from Alexander & Ferguson (1994) For the secondary 
star, we take a mixing-length parameter of a = 2.0 and an effec- 



tive Roche lobe radius given by Eggleton (1983) 



Rl 

a 



0.49^ 



2/3 



0.6q 2 ' 3 +ln(l +q 1 ' 3 )' 



(1) 



where a and q = M2/M1 are the orbital separation and the mass 
ratio of the binary, respectively. 

The mass transfer rate from the donor star via Roche lobe 
overflow can be written as 



2.1. Orbital angular momentum loss mechanisms 
2.1.1. Magnetic braking 

The low-mass main-sequence donor stars in LMXBs were 
braked by the coupling between the magnetic field and the stel- 
lar winds dVerbunt & Zwaan 198 11 1. The loss of specific angu- 
lar momentum in stellar winds is very large because the out- 
flowing ionized particles are bound in the magnetic field lines 
to corotate with the stars out to the magnetospheric radius (see 



Web er & Davis 19671|Me"stel & Spruit l987l|Kalogera 1999| l. In 
LMXBs, the spin angular momentum loss of the donor star 
by magnetic braking causes spin down, while the tidal interac- 
tion between two components would continuously act to spin 
the donor star back up into corotation with the orbital motion. 
Because the spin-up of the donor star consumes the orbital angu- 
lar momentum, magnetic braking indirectly extracts the orbital 
angular momentum of LMXBs. 

A standard magnetic braking model for the angular momen- 
tum loss is widely applied in studying the evolution of cata- 
clysmic variables, and the angular momentum loss rate can be 
written as (Rappaport et al. 1983 ) 



MB 



3.8 x lQ^M 2 Rl{R2lRe) y ^ dyncm 



(3) 



where R2 is the radius, o> the angular velocity of the donor star, 
and 7 is a dimensionless parameter in the range of zero to four. 

However, observations based on rapidly rotating stars with 
a spin period less than two to five days in young open clusters, 
indicate that the standard magnetic braking description should 
overestimate the angular momentum loss rate ( |Queloz 1998 
lAndronov et al. 20031 ). |Sills et al. (2000)| propose the induced 
magnetic braking model, and the angular momentum loss rate 
is given by 



a/2 



J mb -\ \°, ,1/2 



-Kaico 



l\R B M 2 ) 



(4) 



O) > tOr- 



where K = 2.7xl0 47 g cm 2 ( Andro nov et al. 2003). Observations 
show that there is a saturation level in the angular momen- 
tum loss, and w cr j t is the critical angular velocity at which 
the angular momentum loss rate reaches a saturated state. 
Kim & Demarque (1996) suggest that w cr j t is inversely propor- 
tional to the convective turbulent timescale in the star at 200 Myr 
age, 



M 2 = -RMT ■ max[0, (R 2 /R L 



l) 3 ]M e yr-\ 



(2) 



here we take RMT to be 1000 in the calculations. We neglect the 
spin angular momentum of LMXBs because it is smaller com- 
pared with the total orbital angular momentum of the system Q. 
In the calculations we consider two types of mechanisms for or- 
bital angular momentum loss from the binary system, they are 
described as follows. 



If we take typical parameters for the binary, such as Mi = 
1.4 M Q ,M 2 = 1.0 M ,Ri = 10 6 cm,R 2 = 1 R , the spin period of 
the NS is 1 ms, and the orbital period is P m t, = 3 d, the orbital angular 
momentum of the binary is estimated to be J olb = 1.9 x 10 52 gcm 2 s _1 . 
Assuming that the donor star corotates with the orbital rotation, we have 
7, = 7.0 x 10 48 gcmV, J 2 = 27Tk 2 M 2 R 2 2 /P olb = 2.4 x 10 49 gcmV 
(k 2 x OA for main-sequence stars, Dervijo glu et al. 2010} , where J { , J 2 
are the spin angular momentum of the NS, and the donor star, respec- 
tively. It is clear that the orbital angular momentum is three to four 
orders of magnitude greater than the spin angular momentum. 



Went — k>critO" 



(5) 



where w c ri tj0 = 2.9 x 10~ 5 Hz, t , and r are the convective 
turbulent timescales of the Sun and donor star, respectively 
dKrishnamurthi et al. 1997] ) . 

2.1.2. Accretion disk instability and mass loss 

According to the mass of PSR J 17 13+0747 inferred by 
its Shapiro delay, there is a large amount of transfer ma- 
terial that should be lost from this binary system dur- 
ing the mass exchange. Either the accretion disk instabili- 
ties (|Pringle 1981 van Paradijs 1996) or the propeller effect 
( Illarionov & Sunyaev 1975) may be responsible for the mass 
loss even for sub-Eddington accretion. During the mass transfer 
via Roche Lobe overflow, the transferring material would form 
an accretion disk around the neutron star. It is well known that 



Chen & Panei: The progenitor of PSR J1713+0747 



3 



disk accretion in LMXBs may suffer thermal and viscous insta- 
bility, if the effective temperature in the disk is below the hy- 



drogen ionization temperature about 6500 K (van Paradijs 1996 



King et al. 1997 [Lasota 2001). This corresponds to a criti 
cal mass-transfer rate under which the accretion disk is un- 



stable. This critical mass-transfer rate is (van Paradijs 1996 
IDubus et al. 19991 



M„ a 3.2 x 10" 



Mi 



1.4M- 



0.5 



M 2 

I Mr: 



1.4 



M yr 



(6) 



When an accretion-disk instability occurs, the accreting neu- 
tron star will become a transient X-ray source, experiencing 
short-lived outbursts separated by long-term quiescence. When 
the mass transfer rate -M2 is lower than M cr , we assume that 
the neutron star accretes only during outbursts, so the accretion 
rate is M ac f out = -M 2 t r , where f out , t r are the outburst timescale 
and the recurrence time, respectively. Introducing the duty cy- 
cle / = t ou t/t T , we have M ac = -M 2 /f [a typical duty cycle is 
in the range of 0.1-0.001 (King et aL 2003) ]. If -M 2 > M cr , the 
accretion rate is M ac = — M2. In addition, in both cases the mass 
growth rate of the neutron star M\ is limited to the Eddington 
accretion rate Mem, i- e. 



Mi = min[M ac ,M Edd ], 



(7) 



where Medd ~ 1.5 x 10 8 M yr 1 is the Eddington accretion rate 
of the neutron star. 

The transferred matter that exceeds the Eddington accretion 
rate is assumed to be ejected in the vicinity of the neutron star. 
The matter forms an isotropic wind and carries away the spe- 
cific orbital angular momentum of the neutron star. Based on the 
above assumptions, we can derive the angular momentum loss 
rate: 



'iw 



MM 2 
MM, ' 



(8) 



where J = a 2 /iQ is the total orbital angular momentum of the 
LMXB with a circular orbit, fi = MiM 2 /(M l + M 2 ), M = Mi + 
M 2 , and M - Mi -M ac are the reduced mass, the total mass, and 
the mass-loss rate of the binary system, respectively. 

2.2. WD regime: the evolutionary code 

When the LMXB becomes a detached system, the donor star 
evolves into an He-core WD. From this stage, during thermonu- 
clear flashes scenario and in WD regime, we adopt the stellar 
evolutive code LPCODE ( Alth aus et al. 2 005). The code is based 
on a detailed description of the main physical processes in- 
volved in the formation of WD. Radiative opacities are those 
from OPAL (including carbon- and oxygen-rich compositions) 
for arbitrary metallicity (Iglesias & Rogers 1996). These opac- 
ities are calculated for metallicities consistent with the diffu- 
sion predictions. During the WD cooling regime gravitational 
settling leads to metal-depleted outer layers. The equation of 



state is an updated version of that of Magni & Mazzitelli (1979) 



High-density conductive opacities and the various mechanisms 
of neutrino emission are taken from works of |Itoh et al. (1 992) 
Hydrogen burning is taken into account by considering by a 
complete network of thermonuclear reaction rates correspond- 
ing to the proton-proton chain and CNO cycle. Nuclear reactions 
rates are taken from Caughlan & Fowler (1988) We consider the 



evolution of the chemical abundance distribution caused by el- 
ement diffusion during the WD regime. Our treatment of time- 
dependent diffusion is based on the multicomponent gas treat- 
ment presented by Burgers (1969) for element diffusion scheme 
we considered the following nuclear species 'H, 3 He, 4 He, 12 C, 
14 N, and I6 0. Abundance changes are computed according to 
element diffusion and then to nuclear reactions and convective 
mixing. 

In this work, we have computed the complete evolution of 
seven remnants arisen from binary systems as soon as the trans- 
fer of mass has ceased. The masses of these seven pre-WDs are 
listed in Table 1 (M 2 j). The evolution of these WDs is continued 
until they arrive at a luminosity of Log(L/L ) as -4.8. 



3. Results 

To study the characteristics of the progenitor of PSR 
J17 13+0747, we calculated evolutionary sequences of many 
LMXBs. We consider both the angular momentum loss mech- 
anisms described in section 2.2, adopting the induced magnetic 
braking model (eq [4]) for the donor stars with mass less than 
1.5M . Some initial input parameters were set as follows, the 
neutron star has a canonical mass of My = 1.4M Q , the mass of 
the donor star M 2 ^ is in the range of 1.0 - 2.0M o (with a solar 
metallicity of Z=0.02), the duty cycle / = 0.01. 

The characteristic age of PSR J 17 13+0747 estimated by 
Splaver et al. (2005) is t psr ~ 8 Gyr. This implies that the WD 
has a cooling timescale t coo i of approximately 8 Gyr since the 
end of the mass transfer episode. Based on different initial or- 
bital period P 01 -b,i, we stop the binary evolution calculation when 
the binary becomes a detached system and the donor star evolve 
into a WD. Subsequently, a WD evolutionary code (LPCODE) is 
applied to the donor star. If the final donor star mass is in the 
range of 0.25 to 0.3 1M , and the orbital period of the binary is 
approximately of 67.8 d, the corresponding binary was assumed 
to be the progenitor candidate of PSR J1713+0747. 

Firstly, we calculated the evolutionary sequences of an 
LMXB with My = 1.5M Q and P OTb ,i = 2.57 d. In Figs. Q] and 
we plot orbital period and donor star mass, and the mass transfer 
rate and the neutron star mass, respectively. In the early phase, 
the orbital period decrease to be approximately 2.56 d due to 
magnetic braking even if there is no mass exchange. With nu- 
clear evolution, the donor star fills its Roche lobe when its age 
is 2.54 Gyr. When M 2 = O.66M , and P orb = 9.66 d, the bi- 
nary system detaches and mass transfer ceases about 4.9 Myr. 
This detached phase also is found in other studies for the evo- 
lution of LMXBs (Tauris & Savonije 19991 |Han etal. (2000)] 
Podsiadlowski etal. 20021 1. After the temporary shrinkage, the 
donor star fills its Roche lobe again and climbs the giant branch. 
When the age is 2.79 Gyr, the second Roche lobe overflow 
(RLOF) episode terminates. Hereafter, the donor star evolves to 
a helium WD and cools. The evolution end point is a binary pul- 
sar with an orbital period of 68.3 d. The binary consists of an He- 
core WD of O.3O7OM and a neutron star of 1.5184M . Because 
the material is transferred from the more massive donor star to 
the less massive neutron star, the mass transfer proceeds initially 
on short thermal timescale about 10 Myr. In this phase, the mass 
transfer rate is in a range of 10~ 7 to 10~ 8 M o yr _I , and the LMXB 
is a short-lived persistent X-ray source. Subsequently, the mass 
transfer occurs on a much longer nuclear timescale of 0.25 Gyr, 
with a rate of 10~ 9 or lO" 8 M yr"'. This rate is always less than 
the critical mass transfer rate. Therefore, the progenitor of PSR 
J17 13+0747 should be a transient X-ray source about 95% of all 
its life. 
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Fig. 1. Evolution of the orbital period (solid curve) and the donor 
star mass (dashed curve) vs. evolutive timescale for LMXB. 
LMXB has a donor star with an initial mass of 1.5M and an 
initial orbital period of 2.57 d. Horizontal dotted lines on the 
bottom of the figure correspond to the observed mass of WD 
with its error, and horizontal dotted line on the top of the figure 
represents the observed orbital period. 
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Log T eS (K) 
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Fig. 3. Calculated evolutionary tracks in H - R diagram for the 
donor star with an initial mass of 1 .5M Q . The open triangle, solid 
triangle, and open star indicate the beginning of RLOF, the end 
of the last episode of RLOF and the present situation of the donor 
star, respectively. 
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Fig. 2. Evolution of the mass transfer rate (solid curve) and the 
neutron star mass (dashed curve) with the evolutive timescale for 
the same LMXB of Fig. Q] The dotted curve represents the criti- 
cal mass transfer rate, under which the disk accretion is unstable. 
The mass transfer rate has been slightly smoothed for clarity. 



In Fig. [3j we plot the evolutionary path of the donor star in 
the Hertzsprung-Russell diagram (HRD). After the donor star 
evolves into a pre- WD, a hydrogen thermonuclear flash occurs. 
During this stage the element diffusion is a key physical ingre- 
dient in obtaining a thin hydrogen envelope. We stress that ther- 
monuclear flashes also occur in the absence of element diffusion, 
but in this case the residual hydrogen envelope is thick, accord- 
ing to |Benvenuto & De Vito (2004)| and |Panei et al. (2007)1 Fig- 
|4] presents theoretical predicted relation between /° rb,f and Mjj 
for binary radio pulsars. It is clear that our simulation result is 
consistent with the theoretical prediction. 

To investigate the distribution of the initial donor star mass 
and orbital period for the progenitor system of PSR J 1 7 1 3 +0747, 



we calculated the evolution of a large amount of LMXBs for 
different initial input parameters. Table 1 summarizes the calcu- 
lated results. Assuming that the cooling timescale of the WD 
is t coo i = 8 Gyr, as seen in Table 1 makes it clear that the 
donor star with an initial mass of M24 = 1.0 - 1.2M Q would 
require times in excess of the Hubble time in order to evolve 
into a dim WD like the one observed in PSR J1713+0747. 
Therefore, PSR J 17 13+0747 may originated from an LMXB 
with a donor star of M2;, - 1.3 - 1.6M and an initial or- 
bital period of P rb,i = 2.40 - 4.15 d. Thus with an age of 
Tcooi = 8 Gyr, we obtain WD effective temperatures in the range 
of ~ 3870 - 4120 K, approximately in agreement with the ob- 
servational datum 3700+100 K given by Lundgr en et al. (1996)| 
|Benvenuto et al. (2006)| deriv ed a value of 4320 ± 180 K and 
4250 ± 250 K determined by |Schonberner et al. (2000)| In ad- 
dition, the luminosity of the WD was determined to be -4.171 < 
Log(L/L ) < -4.047. Apparently, the final effective tempera- 
ture of the WD in our calculations is sensitive to its cooling 
timescale, and a long cooling timescale can result in a low ef- 
fective temperature. On the basis of these models with real his- 
tory, we estimated the range of surface gravity for the WD in 
7.386 < Log(g) < 7.401 (g in units of g cm" 2 ). Detailed 
tabulations of our WD evolutionary sequences are available at 
http://www.fcaglp.unlp.edu.ar/~panei/papers.html 

Although the final mass of the neutron star (see Nos. 
6 and 7 in Table 1) is slightly higher than the one de- 
it is still in the allowed mass 
by theoretical orbital period - WD mass relation 



rived by Spl aver et al. 2005 
interval 



mass 

(Mi = 1.53^[|ggM ). Though adopting various scenarios, our 
results are approximately consistent with the ones derived by 
|Benvenuto et al. (2006)| 



4. Summary and discussion 

In this paper, we attempt to carefully study the evolution his- 
tory and the progenitor characteristics of PSR J17 13+0747. 
Based on Eggleton's stellar evolution code, we calculated the 
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JNo. 


^2,i 




to 




tf - to 


M 2 , f 




"orb.f 


fp.8 


LogL 8 




L°gg8 




(M a ) 


(days) 


(Gyr) 


(Gyr) 


(Gyr) 


(M Q ) 


(M ) 


(days) 


(Gyr) 


(L ) 


(K) 


(cm s~ 2 ) 


1 


1.0 


5.65 


12.48 


12.68 


0.20 


0.3058 


1.4211 


66.7 


20.68 


-4.0318 


4152 


7.3810 


2 


1.1 


5.09 


8.74 


8.95 


0.21 


0.3069 


1.4220 


68.8 


16.95 


-4.1054 


3999 


7.3909 


3 


1.2 


4.48 


6.12 


6.34 


0.22 


0.3063 


1.4247 


67.6 


14.34 


-4.1349 


3937 


7.3925 


4 


1.3 


4.15 


4.70 


4.92 


0.22 


0.3067 


1.4260 


68.1 


12.92 


-4.1114 


3985 


7.3905 


5 


1.4 


3.70 


3.30 


3.54 


0.24 


0.3068 


1.4964 


68.0 


11.54 


-4.1707 


3872 


7.4001 


6 


1.5 


2.57 


2.54 


2.79 


0.25 


0.3070 


1.5184 


68.3 


10.79 


-4.0471 


4122 


7.3856 


7 


1.6 


2.40 


1.96 


2.20 


0.24 


0.3068 


1.5383 


67.5 


10.20 


-4.1682 


3880 


7.4010 



Note 1. The meaning of the columns are presented as follows. Miy. initial mass of the donor star; P OI b,i' initial 
orbital period; t : age when Roche lobe overflow begins; t f : age when the main Roche lobe overflow ends; t f -t : 
mass transfer timescale; Miy. mass of the WD; Mif. mass of the neutron star; P m \>f. final orbital period; t p $, 
LogL 8 , r e ff 8 and Logg 8 : age, luminosity, effective temperature and gravity when r co = 8 Gyr, respectively. 
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T3 




M 2f (M©) 



Fig. 4. Theoretical predicted relation between P rt>,f an d ^2,f f° r 
binary radio pulsars. The solid curves and dashed curves cor- 
respond to the relations obtained by Rapp aport et al. 1995| and 
Tauris & Savonije 1999| The open star sign represents our sim- 
ulation result. 



evolutionary sequences of a number of LMXBs to fit the ob- 
served data (Mi,M 2 , and P orb ) of PSR J1713+0747. For the 
WD regime we employed the stellar evolution code LPCODE 
dAlthaus et al. 20051 

In calculations, angular momentum loss via magnetic brak- 
ing is included if the mass of the donor star is less than 1.5M . 
The observations show that PSR J 17 13+0747 has only accreted 
less than 0.1 M from the donor star during evolution. Unstable 
mass transfer due to the thermal and viscous instability of the 
accretion disk may be responsible for the extensive mass loss 0. 
Therefore, LMXBs may be transient X-ray sources for most of 
their life. If we adopt a duty cycle of 0.01, the progenitor of PSR 
J17 13+0747 can be constrained to be an LMXB with an initial 
donor star of 1 .3 - 1 .6M (Z=0.02) and an initial orbital period of 
2.40-4.15 d. Adopting a cooling timescale of 8 Gyr for the WD, 
our calculations show that at present its effective temperature is 
between 3870 and 4120 K. Compared with the results derived 
by |Benvenuto e t al. (2006) for Z=0.02, the cooling timescale 

2 Certainly, the propeller effect of the pulsar cannot be excluded 
pilarionov & Sunyaev 1975) 



of the WD in our model is close to the spin-down age of the 
pulsafl We obtained absolute magnitudes of My ^\5. 7+0.4 
and Mi «13.8+0.3 in the V and / bands, and inferred effec- 
tive temperatures T eff * 3301 - 4169 K and T eft - s* 3946 - 4694 
K. Detailed tabulations of our WD evolutionary sequences are 
available at http://www.fcaglp.unlp.edu.ar/~panei/papers.html 

Our accretion disk instability model successfully reproduces 
the orbital period and the mass of the WD. However, the ef- 
fective temperature of the WD is slightly higher than obser- 
vations, although we consider gravitational settling and chem- 
ical and thermal diffusion during the evolution of the WD 
(Panei et al. 2007). Employing Hubble Space Telescope (HST), 



Lundgren et al. (1996) performed optical observations for the 



WD companion of PSR J17 13+0747. Ad opting the tempera- 
ture calibration of (V — I ) colour given by Mo net et al. (1992)] 
they presented the effective temperature of the WD is 
T e s = 3700 ± 100 K. Assuming a blackbody spectrum, 
|Hansen & Phinney (1998)| fitted the HST data, and obtained 
r eff = 3430 ± 270 K. A long cooling timescale of the WD 
may be responsible for the discrepancy between our evolution 
results and the inferred value of observations. In this work we 
take t co = 8 Gyr ( |Splaver et al. 2005) . However, it is very rough 
to estimate the spin-down age of the pulsar from tpsr = P/2P , 
where P and P are the spin period and the period derivative of 
the pulsar, respectively Q We expect that further spectroscopic 
and photometric observations can help us understand the evolu- 
tionary track of PSR J1713+0747. 

As an alternative evolution channel, millisecond pulsars may 
also be formed by accretion-induced collapse (AIC) when the 
mass of the accreting ONe WD exceeds the Chandrasekhar 
mass of 1.4M (INomoto 1984llNomoto 19861 1. Statistical analy- 
sis shows that the birthrate of LMXBs is 1-2 orders of magnitude 
lower than that of millisecond pulsars ( |Cote & Pylyser 1989 
ILorimer 19 95). The AIC channel was raised to interpret this 
birthrate discrepancy ( |Bailyn & Grindlay 1990) >. Based on a sce- 
nario including both the AIC of an accreting WD and unstable 
disk accretion, [Xu & Li (2009)l propose that GRO J1744-28 may 
be evolved from a WD binary. Recently, the estimated birthrates 
of millisecond pulsars show that the often neglected AIC channel 



3 |Benvenuto et al. (2006")l propose that the WD needs 10.8 + 2 Gyr to 
cool when Z=0.02. 

4 The standard spin-down age of the pulsar is r PSR = ^^ [1 - 

(Po/P)" where P is the initial period and n the braking index of 
the pulsar. When Po <k P and n = 3, tpsr = P/2P. However, the de- 
rived braking indice of several young radio pulsars are all less than 3 
\Xu & Qiao 200T1lChen & Li 2006b . 
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cannot be ignored ( Hurley et al. 20 IP") . Therefore, the formation 
possibility of PSR J17 13+0747 via the AIC of a massive WD 
cannot be ruled out. We plan to explore the AIC channel for the 
formation and evolution of PSR J1713+0747. 

Of course, there are many uncertainties in our input physics. 
Firstly, the orbital angular momentum loss plays a vital role 
in the evolution of LMXBs. Some mechanism of extracting 
angular momentum possibly exists besides magnetic brak- 
ing and mass loss due to the accretion disk instability. For 
example, the propeller effect of the neutron star may carry 
away the specific orbital angular momentum at the magne- 
tosphere of the neutron star ( jlllariono v & Suny aev 1975| l. 
In addition, a circumbinary disk around LMXBs may also 
extract orbital angular momentum efficiently from the bi- 
nary system ( |Spruit & Taam 200~T] Taam & Spruit 2001} . 



Secondly, our results strongly depend on duty cycles, which 
are poorly known. In this work, we take a constant duty 
cycle. However, duty cycles may change from one system to 
another (Lasota 2001 ), and change with the orbital period and 
the mass transfer rate. Finally, a gamma-ray flux caused by 
the interaction of the magnetic field of the neutron star with 
the accretion disk may drive a wind from the donor star by 
heating its outer layers; i. e., the donor star may be evaporated 
(|van den Heuvel & van Paradijs 1988] IKluziniak et al. 19881 
IRuderman et al. 1989a||Ruderman et al. 1989bl >. 
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